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I. INTRODUCTION 



One-dimensional lattice models with incommensucate hopping elements or on-site potentials have been 
studied for many years from different points of viewld'EI. Many unusual properties of the quantum spectra 
and wave functions have been discovered for various kinds of aperiodicityEl. Physically, such modelsjhave 
applications to several problems such as incommensurate crystajsQ, semiconductor heterojunctionsQ, the 
incommensurate phase of spin-Peierls systems such as CuGeOaQ, and the Azbel-Hofstadter problem of 
particles hopping on a two-dimensional lattice in the presence of a magnetic fieldQ. If the incommensura- 
tion is close to dimerization (wave number tt) , then the models have an additional interest in thepcontext 
of metal-insulator transitions near half-filling and spin-Peierls systems near zero magnetizationEliJ. 

In this paper, we study a model of fermions with a hopping which has a weak sinusoidal incommensurate 
term close to dimerization. We will show both numerically and analytically that this has some unusual 
properties, particularly inside the gap which exists exactly at dimerization. The major surprise is that 
states appear inside the gap as soon as we move away from dimerization. These states carry a finite 
weight, and therefore contribute to quantities like the specific heat at low temperatures, as we will show. 
In Sees. II A and II B, we discuss free fermions or the XY spin-1/2 chain close to dimerization (and 
half-filling) for which quantitatively accurate results can be obtained by analytical methods. We will use 
a continuum Dirac field formulation, followed by perturbation theory and a WKB analysis, to find the 
locations and widths of the energy bands inside the q = gap; here q is the deviation from tt of the 
wave number of the incommensuration. In Sec. II C, we will explain how the limit g — > get reconciled 
with the results for g = in a somewhat different model, namely, the Azbel-Hofstadter problem in 
which electrons hop on a two-dimensional lattice in the presence of a magnetic field. In Sec. Ill, we use 
bosonization to argue that similar differences between g — > and q = occur for interacting fermions or 
the XX Z spin-1/2 chain, which is the more general and interesting case. In Sec. IV, we show that such 
peculiarities do not occur in the vicinity of any wave number which is different from tt. In Sec. V, we 
consider continuum theories which go up to higher orders in the strength of the incommensuration; this 
enables us to discuss the complete energy spectrum of fermions with a weak incommensurate hopping. 
In Sec. VI, we make some concluding remarks. 

Our work differs from earlier ones which have concentrated on the effects of infiemHiensuration which 
are either large or close to "highly" irrational numbers such as the golden ratioEfHa^EI. Further, we will 
use techniques from continuum theory which have not been used much in this area (except for the recent 
work in Ref. 8). 



II. NONINTERACTING FERMIONS CLOSE TO DIMERIZATION 

A. Numerical results 

We begin with the following Hamiltonian for noninteracting and spinless fermions on a lattice 

H ^ - Jn C„+i + ct+iC„ ) - ^1 ^ c]^ Cr, , 

n n 

J„ = 1 + (5 cos(7r -I- q)n , (1) 

where we will assume that 6^1 and q ^ tt. (In that limit, we will show later that it does not matter if 
the incommensurate term is put in the chemical potential fi (i.e., on-site) rather than in the hopping J„). 
We set /i = as we are interested in energies close to zero. If 5 = 0, the model can be easily solved; the 
dispersion relation is E{k) — — cosfc. In the ground state, all the states with momenta lying in the range 
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[—kp,kp] are filled, where the Fermi momentum kp = tt/2. The Fermi velocity is equal to sinkp = 1. 
(We will set both the lattice spacing and Planck's constant equal to 1). 

If 5 ^ but q = 0, i.e., with dimerization, the model can still be solved analytically. There is an energy 
gap extending from —S to 6. Let us now consider nonzero values of q. [Since cos{Tr+q)n = cos{qn) cos(7rn), 
we see that a small value of q is equivalent to dimerization term whose amplitude S cos{qn) varies slowly 
with n]. For any rational value of q/i: — AI/N, where M and N are relatively prime integers, we have a 
periodic system with period P equal to N ii M + N is even and 2N is M + iV is odd. The one-particle 
spectrum of Eq. (Q) can be found by solving the discrete Schrodinger equation 

[Jn V'n+l + Jn-l V'n-1 ] ^ E ^pn , (2) 

which is equivalent to the equation 

/ Tpn+l \ _ { -2E/Jn -Jn-l/Jn \ { i>n \ /q\ 
Un J " V 1 J { J ■ 

The allowed values of energy can therefore be found by computing the transfer matrix M{q, E) obtained 
by multiplying together P matriceaj, 

M{q,E) = n {'^^/'^" -^-i/^") . (4) 

n=l 

Since detM{q, E) — 1 and txM{q,E) is real, it is clear that the eigenvalues of M{q,E) must be either 
of the form r and 1/r, or of the form exp(ir) and exp(— ir), where r is real in either case. In the 
former case, the wave functions ipn diverge exponentially as n — > cxd or —oo; hence they are not allowed 
physically. In the latter case, the wave functions do not diverge at infinity and are allowed. Thus, if 
|tr M{q, E)\ = 2\ cosrj < 2, the energy E is allowed in the spectrum; otherwise that energy is not allowed. 
By following this method and sweeping through a large number of values of q and E, we obtain the picture 
of energy bands and gaps (shaded and unshaded regions, respectively) shown in Figs. 1 and 2, taking 
5 = 0.05. [We have sampled the energy values with a resolution of dE — 10~^. It may be useful to 
clarify here that the numerous vertical gaps in Figs. 1, 2, 5, 6 and 7 have no physical significance; they 
are present because we have omitted the regions ot q/n whose rational approximations have such long 
periods P that a numerical computation would have taken too long]. We have scaled both E and q in 
units of S because the pictures turn out to depend only on those two ratios. We immediately see that 
the pictures look much more complicated than the situation in which q is exactly equal to zero; in that 
case, there is precisely one gap extending from E/S = to 1. [We will consider only positive values of q 
and E since the spectrum is invariant under either q ^ —q (as is clear from (^) ov E —E; the latter 
is clear from the invariance of \tTM{q, E)\ under an unitary transformation by the Pauli matrix 

'' = (J-i) (^) 

which reverses the sign and simultaneously converts E — > —E for each of the P matrices appearing in 
(^). Thus the energy bands shown in Figs. 1 and 2 should be understood as continuing to negative values 
of both E and q by refiection about those two axis] . 

Fig. 1 shows that as q/6 increases, the gaps shrink rapidly. We will show that this can be understood 
using perturbation theory to n^^ order, where n is an odd integer, and that the gaps shrink as (S/q)". We 
have therefore labeled the gaps in Fig. 1 by the integers n = l,3,5,.... More interestingly, we observe that 
all the bands approach the origin (q, E) = (0, 0). The widths of the low-lying bands vanish exponentially 
fast as q/6 ^ as we will argue from a WKB analysis. This makes it impossible to find these portions of 
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the bands by looking for energies satisfying |tr A'I{q,E)\ < 2 using any reasonable energy resolution dE. 
We therefore find these thin portions of the bands by looking for minima of trM'^{q, E) as functions of E. 
If the energy resolution was infinitesimal, these minima would be found at trM^(q, E) = 2cos(2r) = —2, 
i.e., trM{q,E) — 2cosr = 0, which occurs inside an energy band. In practice, due to the finite energy 
resolution dE, these minima numerically yield single points in energy which lie within a distance of dE 
of a band. We are thus able to find the locations of the thin regions without having to use a resolution 
smaller than dE = 10~^. In Fig. 2, we show these points (thin regions) which smoothly join on to the 
wider regions of the bands. We also show six curves which are the analytical results of a low-energy 
theory of the model which involves solving a 1 + 1-dimensional Dirac equation in a periodic potential. We 
will also prove that the number of bands in the region < E/6 < 1 increases a.s 26 /nq a.s q/S ^ 0, while 
the number of states in each band is equal to q/n when normalized appropriately in the thermodynamic 
limit, i.e., the number of sites L oo. This will show that the number of states lying inside the q = 
gap is finite and equal to 26 /tt'^ in the limit q ^ 0; this implies that g = is a rather singular point. 



B. Analytical results 

To begin, let us write the Hamiltonian in (Q) as a sum H = Hq + Vs, where Vs is the incommensurate 
term 

Vs = - 2 XI cos(gn) ( cicn+i + ) , (6) 

n 

where Q — tt + q in this section. If 



|fc) = —= exp{ikn) \n) (7) 



denotes a momentum eigenstate of Hq, then the matrix elements of Vs are given by 

{h\Vs\h) = - ^ ( SQ+k,M + h+k.M ) ( e-^'^ + e^'^^ ) ■ (8) 



Let us understand the gaps for large values of q/6 by using perturbation theory about 6^0. If q 
and 6 are both small, the states close to zero energy are dominated by momenta near ±7r/2. Since 
the incommensuration term Vs in has Fourier components at momenta ±(7r + q), we see that the 
gaps above zero energy result from the breaking of the energy degeneracy between the two states with 
momenta ki = 7r/2 + ng/2 and k„+i = —Tr/2 — nq/2, where n = 1,3,5, .... These two states are connected 
to each other through the n ~ I successive intermediate states with momenta /c2 = — 7r/2 + (n — 2)q/2, 
fca = tt/2 + (n — 4)q/2, ^4 ~ — 7r/2 + {n — 6)q/2, fc„ — 7t/2 — [n — 2)q/2, since we must have 
/cj-fi — kj = ±(7r + q) mod 27r. Since all the matrix elements {kj^i\Vs\kj) are approximately equal to 
±.16/2 (for j ~ 1, 2, n) for small values of q in Eq. (^), and all the energy denominators Ej — E\ are of 
the order of q (for j = 2, 3, n), we see that the two-fold degeneracy (at the energy E\ — En+i — nq/2) 
gets broken at the n**^ order in perturbation theory to produce a gap of the order of J"/?"^^- To be 
explicit, we find that the upper and lower edges of the two gaps labeled by n = 1 and 3 in Fig. 1 are 
given by 

E±{1) ^ ^ , 1 

5 26 2' 

6 26 ^ 16q 32g2 ' 
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and we have verified that this agrees weU with the numericaUy obtained boundaries. We can show that 
the general formula for the gap AE{n) — E^{n) — E_{n) labeled by the odd integer n = 2p + 1 is given 

by 



to lowest order in 5. 

We can count the number of states in the band lying between gaps n and n + 2 as follows. For 
(5 = 0, let us normalize the number of states so that the total number is 1; since the momentum k goes 
from — TT to TT in that case, the density of states in momentum space is l/(27r). If we now turn on a 
very small value of 5, we find that the band lying between the gaps n and n + 2 is made up of linear 
combinations of the states lying between the two momenta intervals [7r/2 + nq/2, -k 12 + [n + 2)g/2] and 
[— 7r/2 — ng/2, — 7r/2 — (n + 2)(7/2]. The total number of states in these two intervals is 2((7/27r) = qj-K. 
Thus each band contains qjix states. Now, this number cannot change if we change 5, and we can therefore 
use the same number below in the opposite limit where qjS \s small. 

Let us examine the more interesting regime where qjb is small and Ejb < 1. For analyzing this, it is 
useful to consider the continuum limit. In this limit, the Fermi field tpi^) — Cn can be written as 

ip{x) = i>R{x) exp(^) + i}L{x) exp(-^) , (11) 



where V'fl a-nd '^h denote the right- and left-moving fields, respectively; they vary slowly on the scale of 
the lattice spacing. We substitute ( pi] ) in Eq. (|l]) and drop terms like exp(±i7ra;) which vary rapidly. We 
then find the following Dirac-like equations for the two time-dependent (Heisenberg) fields 



i ( djdt + djdx ) ^pji — iS cos(qx) tAl = , 
i ( d/dt — d/dx ) ifjL + iS cos{qx) tpR ^ . 



(12) 



The energy spectrum can be found from (112) by defining ^± = ^pn ± ipL which satisfy the equations 



[ - d^/dx^ + COS^iqx) - Sq sm{qx) ] V+ , 
[ - d^/dx^ + cos2(ga;) + 6q siii{qx) ] V- = V'- ■ 

It is sufficient to solve one of these equations, say, for 7/;+, since tp^ is related to tp+ by 



(13) 



^-{x) = — —d/dx + S cos{qx) ip+ix) , 
El i 

provided E ^ 0. The energy spectrum can be found by solving the time-independent equation 

[ - d^/dx^ + cos^iqx) - Sq siniqx) ] ^1;+ = E^ V+ . 



(14) 



(15) 



Eq. (|T^) has the form of a Schrodinger equation (with "energy" E'^) ia the presence of a periodic 
potential. The potential is similar but not identical to Mathieu's equationliJ. We therefore again expect 
bands to form. To simplify the notation, let us shift x by TT/2q and then scale it by a factor of q to make 
the period equal to 27r. We then get 

[ — d^/dx^ + sin^ X — - cos a; 1 ^+ = ■0+ • (16) 

q 
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By Floquet's theorem0, the solutions must satisfy il^+{x + 2ti) = e''^ip+{x), where 9 goes from to 
±7r from the bottom of a band to the top. We observe that there is an exact zero energy state with 
i'+ix) = exp( I cos a; ) and ^-(x) = for which 9 = 0. In general, the nonzero energy states can only be 

found numerically. But if q/S is large, the positions of the low-lying bands (with E'^/q'^ <C S/q) can be 
found analytically as follows. For energies much lower than the maxima of the potential, we can ignore 
tunneling between the different wells to begin with. We note that there are two kinds of wells in Eq. (|T^); 
those centered about x = 2m7r and those centered about x = (2m + l)7r, where m can be any integer; the 
former wells are deeper than the latter. However, the perturbative energy levels (i.e., ignoring tunneling) 
are identical in the two types of wells, except for the E = state which only exists in the deeper wells. 
This identity follows from Eqs. (1^) and ( p^ which show that corresponding to any (nonzero) energy 
eigenstate ip+i^) which lies in a deeper well, there is an eigenstate ip-{x) which lies in a shallower well. 
The n"' energy level in a deeper well has the same value as the (n — l)*'^reiiergy level in a shallower well, 
for n > 1. (This is an example of supersymmetric quantum mechanicai^) . It is therefore sufficient to 
find the energy levels of (^6|) which lie in one of the deeper wells, say, at a; = 0. Near the bottom of that 
well, we have a simple harmonic potential with small anharmonic corrections. Ignoring the anharmonic 
terms, we find the energy levels to be simply given by 

El = 2nqS , (17) 

where n = 0, 1, 2, .... Next, we include the anharmonic corrections perturbatively. Up to second-order in 
perturbation theory, we get the more accurate expression 

En = J2nqb 1 — . 18) 

^^^86 128 (52 J ^ ^ 

In Fig. 2, we show the six curves corresponding to n = 1 to 6 in Eq. (p^); for n = 0, we simply get 
a straight line lying at zero energy. We see that all the curves agree extremely well with the numerical 
data, even up to i?„ of the order of d where the harmonic approximation breaks down and the band 
widths become noticeable. 

We will now consider tunneling between wells using the WKB method. For n <^ 6/q, the energy levels 
E^/q^ of (^6|) lie in small intervals centered about x = rmr, where m can be any integer for n > I but 
must be an even integer for the n = band. In these small intervals, the potential is simple harmonic to 
a good approximation, with the harmonic frequency being equal to 

b=S/q. (19) 

We will assume that b is large so that the band widths are small. The splitting of the n^^ band is given 
by the expressionllj 

A{El) 4b 



- exp( - / dx y/V(x) — 2nb ) 

u2 „;„2 



V{x) = ¥ sin-" X - b cosx , (20) 

where xi and X2 are the turning points, and we have approximated E^/q^ by the lowest order result 2nb 
given in (|7|). For n> 1, we take the turning points to be Xi ~ ^ {2n + l)/6 and X2 — tt — \J {2n — l)/6, 
while for n = 0, xi ~ 1/Vb and X2 — 2t: — l/\/b. To lowest approximation, the exponential factor in ( |20| ) 
is given by exp(— 26) for n > 1 and by exp(— 46) for n = 0. However, since A(£'2) ~ 2£^„Ai?„ for a small 
width AEn if n > 1, but A{E^) = (AEq)^ (since Eo = 0), we see that the splitting AE^/q ~ exp(-26) 
for all n. We can do a more accurate calculation to find the prefactors; we find that 

AEn 23"+i 2S 

exp( ) (21) 

.gr/ q 
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for all values of n. We should note that the expression in ( pT| ) for n — Q actually gives half the width of 
that band since the n = band extends an equal amount to negative values of energy. 

Eq. (|2^) explains why, with a numerical method using a finite energy resolution dE, the bands rapidly 
become thin and shrink to isolated points as g — s- in Fig. 2. In Fig. 3, we compare the numeri- 
cally obtained widths of the lowest four bands (n = 0,1,2,3) with the WKB expressions in (pT|). The 
points in that figure indicate the thirty smallest values of width that we could find using the condition 
|tr M{q,E)\ < 2 with the energy resolution dE = 10^^; the solid lines indicate the WKB expressions. 
We see that the agreement between the two becomes worse as q/d increases, particularly for the larger 
values of n. This is expected because the WKB method is accurate only if n <C 6/q. 

We will now prove that the number of states lying below the line E/6 ~ 1 is proportional to 6. For 
large S/q, we can use a semiclassical phase space argument to count how many states lie below any 
given energy. The "Hamiltonian" on the left hand side of Eq. ( [T5| ) has the form + V{x), where 
V{x) — (5^ cos^(qx) — Sqsin{qx). Hence the number of states up to energy E is given by the phase space 
integral 

Iff dxdp ^, 2 2 



""^^^ ^ L J J ^^^^ ^^"^^ ^ ' ^^^^ 

where we have divided by the length L of the system for normalization; the 0-function is defined to be 
1 and if its argument is positive and negative, respectively. On doing the momentum integral in (|2^), 
we get 

i^iE) = ^ j^j'' dx ^E^ ~ V{x) e( E^ - V{x) ) , (23) 

where we have used the fact that V{x) has period 2-K/q. For small values of q/5, we then see that the 
number of states lying between E = and 5 is 

^-5) = % , (24) 



which is independent of q. Incidentally, this also yields the number of bands lying below E — 5. Since 
each band contains q/ir states, Eq. ( p4| ) gives the number of bands to be 25 /-Kq. We have verified that 
this estimate agrees very well with our numerical results. 



We can obtain the density of states p{E) from the expression in (|23|) by differentiating it with respect 
to E. In the limit q/5 ^ we get 



P{E) ^ - drj 



1 - (|)%inS^ 



-1/2 



where U = sin"^(S/(5) if E<5 and 7r/2 if E>S 



(25) 



These are shown by the solid lines in Fig. 4; note the linear behavior near E/S — and the logarithmic 
divergence as i? — > (5 from either side. This is to be contrasted with the density of states found exactly at 
g = 0. In that case, p{E) vanishes for E < S and is equal to l/[7r-y/l — (S/E)^] for E > 5. This is shown 
by the dashed line in Fig. 4. 

The difference in the spectrum for q = and g — + should show up most clearly in the specific heat at 
temperatures T <^ 5. Since the chemical potential is zero, the free energy per site is given by 

F = -2T dE p{E) ln( 1 + e"-^/"^ ) , (26) 

"'0 
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The factor of 2 on the right hand side is because we have to sum over both fermions and holes which 
have the same density of states. (We set the Bohzmann constant equal to 1 for convenience). Then the 
specific heat per site is 

Cv = -Ty. (27) 



For q — 0, the energy gap from —S to +6 means that the specific heat vanishes exponentially at low 
temperature. But for q — > 0, Eq. ( p5| ) shows that the density of states goes linearly as p{E) ~ E/{ii5) 
for E <^ S. This implies that the low-temperature specific heat goes as T'^/S which is very different from 
the exponential behavior at q = 0. 

We would like to point out here that the above results remain unchanged if the incommensurate term 
is present in the chemical potential fi rather than in the hopping J„ in Eq. (^, provided that S and q 
are both small. In other words, if we start with the Hamiltonian 

H = - ^ 51 ( ^" "^"+1 + cl+ic„ ) - 6 cos(7r + q)n ^ 4 c„ , (28) 

n n 

we get the same continuum theory as discussed above, if we define the Fermi fields as in Eq. (|l^), followed 
by a phase redefinition iJjl i'tpL- This equivalence between the two continuum theories will be used 
below. 



C. Reconciling q — with q ^ 0: the Azbel-Hofstadter problem 



Let us ask: why does the model in Eq. (|lj) show such different behaviors in the limit q — and 
at (7 = 0? To answer this question, it is useful to generalize the form of the incommensurate hopping 
in (|^) from (5cos(7r + q)n to Jcos[(7r + q)n + 77], where 77 is a phase lying between and 2tt. If q/ir is 
irrational (this is the generic case to consider if q is nonzero), the spectrum does not depend on rj. [If 
q/n = M/N is rational, the spectrum does depend on 77, but it varies less and less with 77 as — > 00. 
In our numerical computations, we introduced a random phase rj for each value of q and found that this 
had no noticeable effect on the two figures] . This insensitivity to 77 can also be seen from the continuum 
theory in Eq. ( [l^ where a potential of the form cos{qx + 7/) can be transformed to cos{qx) by shifting x 
appropriately. However, exactly at g = 0, the hopping has a term like (5cos77COs(7r7i) and the spectrum 
depends significantly on 7/; for instance, there is a gap from = up to |(5cos77|, and the number of 
states up to an energy E is given by 

v{E,ri) = - ^JE"^ - d^cos'^T] (29) 

TT 



for |(5cos77| < E ^ 1. We may now consider taking an average oi this number over all possible values of 
77. Thus, the number of states up energy E is given by Eq. (|29| ) to be 

v{E) = ^ ^/E-^ - 52 cos2 7; e( ^2 _ ^2 ^Qg2 ^ ) (3Q^ 

Jo 27r2 

This agrees with (^3|) if we substitute rj = qx. Thus, the limit g — > agrees with the point q = Q provided 
that we average over 77 in the latter case. 

The question now arises: is there a physical system where an average over the phase 77 occurs naturally? 
Some experimental systems sitting at g = (the dimerized point) are more likely to choose a particular 
value of 5 with 77 = 0, rather than average over many values of 77; the limit (7 — > and q exactly equal to 
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zero will therefore differ from each other. For simplicity, we will choose rj = in the following sections to 
describe such systems. 

There is, however, a problem in two dimensions which averages over rj in a natural way. This is the 
Azbel-Hofstadter (AH) problem in which electrons hop on a regular lattice in the presence of an uniform 
magnetic field pointing normal to the plane. Let us consider a square lattice with the sites being labeled 
by two integers (m,n). (We set the lattice spacing equal to 1 as usual). The hopping element from the 
site (m, n) to the site (to+ 1, n) will be denoted by t(rn, n; x), while hopping from (m, n) to (m, n+1) will 
be denoted by t{m,n]y); these will be complex numbers in general. Suppose that the hopping elements 
are all real in the absence of a magnetic field, with i(m, n; x) = ti and t(m, n; y) = ^2- In the presence of 
a magnetic field, let $ denote the magnetic flux through each square. Then some of the hopping elements 
must necessarily be taken to be complex, with the sum of the phases of the four hopping elements taken 
anticlockwise around each square being equal to Q = e^/{hc), where e is the electronic charge and c 
is the speed of light. We can interpret Q/{2tt) as the ratio of the flux through each square to the flux 
quantum $0 = 2T:hc/e. Clearly, Q and Q + 27r are physically indistinguishable; thus, Q = tt corresponds 
to the largest magnetic field that can appear on a lattice. Now, there are many choices possible for 
distributing the phases amongst the different hopping elements; the different choices are related to each 
other by gauge transformations. Let us choose the Landau gauge in which the hopping elements in the 
y direction stay real and are equal to t2, while the hopping elements in the x direction are given by 

t(m, n; i) = ti exp(— iQn) . (31) 

The lattice Schrodinger equation satisfied by the single-particle wave functions is then given by 
ti ( exp{iQn) 4)m+i,n + exp(-i(5n) 

m— l.n ) + ^2 ( (32) 

Since the hopping elements do not depend on the ai-coordinate m, we may assume the wave functions to 
have the factorized form 

0m, « = exp(iTO77) ipn , (33) 

where the momentum rj can take any value from to 2tt. Eq. (|3^ ) then takes the form 

1 ti e 

- o ( V-'ri+i + V'n-i ) - — cos(Qn + 7;) = _ ^„ . (34) 

Z 12 zt2 



This is the well-studied Harper's equation with an on-site incommensurate term. Further, we see that a 
phase r] appears naturally in that term, and that we have to consider all values of rj lying in the interval 
[0,27r]. 

Let us consider the extremely anisotropic limit S = ti/t2 << 1- We will show that this problem has 
the same continuum limit as the problem in which the incommensurate term is in the hopping. If we 
define 

S ^ — and E ^ — , (35) 

t2 2t2 ^ ^ 

then the Hamiltonian corresponding to Eq. ( |3^ ) can be written as H = Hq + Vs , where the incommen- 
surate term is 

Vs = - (5 ^ cos{Qn + rj) cj,c„ . (36) 



The matrix elements of Vs between momentum eigenstates of Hq are given by 
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(37) 



For Q close to tt and fci, k2 close to ±7r/2, this has, up to factors of i, the same form as (g), except that 
we have to average over 77 in the AH problem. We thus see that the anisotropic limit of the AH problem 
provides us with a physical realization of Eq. (|l|) except that we have to consider a phase average of the 
incommensurate term. Dimerization in the one-dimensional problem corresponds to Q = tt, namely, the 
magnetic flux through each square of the two-dimensional lattice is equal to half a flux quantum. We 
may therefore carry over all our results to that particular limit of the AH problem. For instance, the 
density of states for Q ^ tt has the form given by the solid lines in Fig. 4. Note that the energy spectrum 
does not depend on the phase 77 for generic (i.e., irrational) values of Q in Eq. (^^. So if denotes 
the number of lattice sites in the x-direction, the number of values of rj is also equal to since 77 takes 
all values from to 2tt in steps of 2tt/Lx- Hence the density of states for the AH problem is equal to 
times the density given in Fig. 4 for the one-dimensional problem. 

We note that the energy eigenvalues in Eq. (|3^) are invariant under the exchange ti <-^- ^2 corresponding 
to an interchange of the x and y hopping elements in the Landau gauge. Thus the AH problem has the 
duality symmetry S — > 1/S. Our problem with an incommensurate hopping has no such symmetry. 

With the present day expertise in fabricating quantum dot arrays, it may be possible to construct a 
two-dimensional lattice whose spacing is large enough that the flux through each square could be made 
equal to half a flux cmantum with currently available magnetic field strengths. (A similar suggestion was 
made by Hofstadteitl) . Consequently, our expression for the density of states of such a system could be 
experimentally tested. 



III. INTERACTING FERMIONS CLOSE TO DIMERIZATION 



Let us now return to one dimension. In this section, we will argue that the difference between small q 
and q — (in particular, the presence of states within the q — gap) persists for the more interesting case 
of interacting fermions, i.e., for Luttinger liquids. To this end, consider adding a four-fermion interaction 
term like '^"'^"'-I+i'^n-i-i to the Hamiltonian (Q). Equivalently, we can consider an XXZ spin-1/2 
chain governed by 



[ ( 1 + Scos{n + q)n ) ( S^S^^, + Sy^+i ) + D S^S^^, ] . (38) 



which is related to the interacting fermion theory by a Jordan- Wigner transformationt^l. If the incom- 
mensurate |tarm is absent (S = 0], it is known from the exact Bethe ansatz solution and conformal 
field theorjiia that the model in (Bq) is gapless with a linear dispersion of the low-energy excitations if 
— 1 < D < 1. If we now add a dimerization {q — and S is small), then the system becomes gapped and 
the gap scales as 

AE = , 

where K = ^- , (39) 

2cos-i(-i:') ' ^ ^ 

where < cos~^(— I?) < tt. [Actually, a gap opens up only if K < 2. If K > 2, the dimerization term 
is irrelevant in the sense of the renormalization group, and a gap is not generated]. It is useful to state 
the result in Eq. ( |39| ) in the language of bosonization. We introduce a bofionic field 4>{x, t) such that 
bilinears in Fermi fields have local expressions in terms of the bosonic fieldlla. For instance, 

V'lji/'L ~ exp(i20) , 

- exp(-i20) , (40) 
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and 



d(j)/dx 



TT [ p{x) - po 



(41) 



where p{x) is the local fermion density, and po is the average density. For q = 0, the model in (38) is 
equivalent, at low energies, to a sine-Gordon theory described by the Lagrangian densityllj 



1 

2'KvK 



{d(f>/dtf - {dc^/dxf 



(42) 



where v is the velocity of low-energy excitations and a is a positive constant; their numerical values are 
not important here. The main point is to note the exponent of 5 in the coefficient of,cos(2(/)). (There are 
also factors of In 5 due to the presence of a marginal operator in the XX Z modeO, but we will ignore 
such terms here). The number K and the velocity v are the two important parameters which determine 
the low-energy, long wavelength behavior of a Luttinger liquid. 

It is important to observe that the incommensurate term which is linear in 5 in the original fermionic 
theory in (|38| ) becomes a cosine term with a different exponent for i5 in the low-energy bosonic theory. 
This is because a nontrivial renormalization xiccurs in the process of deriving the low-energy bosonic 
theory from the microscopic fermionic theoryuilj. This renormalization occurs even if the fermions are 
noninteracting, i.e., for the XY spin-1/2 chain with D = in ( ^Sf ) and K — \; the reason for this is that 
the sine-Gordon theory is always strongly interacting. These strong interactions are also responsible for 
the large renormalizations of the correct quantum spectrum compared to the naive (i.e. classical) spec- 
trum that one obtains from the sine-Gordon theory, namely, the classical soliton mass for thp fermionic 
excitations and the quadratic fluctuations around spectrum for the bosonic excitationtO. The non- 
interacting fermionic model in (|^) has no such renormalizations, which is why we did not use bosonization 
in the earlier part of this paper. 

Let us now change q slightly away from zero. Since cos(7r -I- q)n = cos(7rn) cos(qn) on a lattice, we 
have a dimerization whose coefficient 5 cos qn varies very slowly over the system. Thus, the low-energy 
excitations, which had a relativistic dispersion with the mass now have a space- 

dependent mass \5 coa{qx)\^/'^'^~^\ To put it differently, bosonization yields a theory of the sine-Gordon 
type, except that the coefficient of the cos(2(/)) term in Eq. ( |4^ ) gets modified to produce 

- a |(5cos(ga;)|2/(2--ff) [ i _ cos(2(?!)) ] , (43) 

Unlike ([42|), the theory in Eq. ( [43| ) cannot be solved analytically, either in classical or in quantum 
mechanics. However we can make some qualitative statements about the low-energy spectrum ii q/5 <^ 1. 
When calculating the spectrum of small oscillations about a classical ground state (p = 0, we find it 
convenient to first shift x by 7r/(2g) to change cos(gx) to sin(ga;), and to then scale x by a factor of 

a = (g<5)i/(3--^) . (44) 



C = 



i 

IttvK 



{d(j)/dtf - {d<i>ldx) 



This gives the eigenvalue equation 

_ „2 ^ + ^^Kav = {-^ <l> (45) 

ox^ \ a ) 

for eigenstates lying close to the origin a; = 0; we have approximated sin(qx) by qx and sin(20) by 2(/). 
Since ( ^5| ) is the Schrodinger equation with a confining potential, we see that the energy can take several 
discrete values which are given by numerical factors multiplying a. [Note that our earlier results for 
noninteracting fermions agree with this scaling argument if we set K = \\. These discrete values will 
then spread out into bands once we include tunneling between the different wells centered at the points 
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X — mn. Similarly, we can find the energy of a classical soliton which goes from = at x — > — cxd to 
(j) = IT at X oo; once again we can show by scaling that this will be given by a times some numerical 
factor. [Of course, all the numerical factors will get renormalized due to quantum corrections, but the 
power-law dependence on q6 is not expected to change] . We therefore see that all the low-lying excitations 
have energies of the order of {qSy/^^~-^'>; this is much smaller than the gap which exists exactly at g = 
given by Eq. (|3^), since we are assuming that q/S is small. We thus see that in the limit 5 — > 0, there are 
many states which lie within that gap. Indeed, we can use WKB quantization for Eq. (^) to estimate 
that the number of bands lying between zero energy and the gap AE in ( ^9|) is of the order of AE/q, 
in the limit q ^ 0. The number of states in the same interval is therefore of the order of AE; this also 
follows from the semiclassical argument presented below. 

The semiclassical phase space estimate given in Eq. can be used to find the density of states p{E) 
at energies much smaller than the q = gap. If i5 <C S^^^-^\ the wave function for that state lies in a 
region where we can approximate shi(qx) by qx. Ignoring various numerical factors, we then get 

p{E) = 4- — \Ie^ - {Sqx)y<-^-^) e( E^ - (,5<zx)2/(2-^) ) ^ . (46) 

dE Jq t: ^ d 

Since the chemical potential for these bosonic excitations is zero, the free energy per site is given by 

/>CXD 

F = T dE p{E) ln( 1 - e"^/^ ) . (47) 

Jo 

The specific heat at temperatures much lower than the q = gap therefore scales as T^~^ which is 
again much larger than the exponential dependence which occurs at g = 0. [It is tempting to compare 
this exponential versus power-law dependence of the specific heat with the forms of the specific heats 
of the niagnetic excitations observed in the dimerized {q = 0) and incommensurate (q ^ 0) phases of 
CuGeO^. However that system is much more complicated because it has strong phonon-spin couplings; 
we have ignored such terms in our model by assuming that the incommensuration is static] . 

We should note that Eq. (|2^) counts the number of one-particle fermionic states, and the specific 
heat calculated from that gets a contribution from states with all possible fermion numbers. Eq. (|4^), 
however, only counts the number of one-particle bosonic states, i.e., fermion-hole excitations. Hence 
the specific heat calculated from that only gets a contribution from states with zero fermion number. 
However the complete expression for the specific heat is expected to differ only by numerical factors from 
the bosonic one, and should therefore have the same kind of exponential or power-law dependence on the 
temperature. 



IV. FERMIONS NEAR OTHER COMMENSURATE FILLINGS 



It is interesting to note that the peculiar difference between the limit q and 5 = only occurs 
near dimerization, i.e., near wave number tt. There is no such singularity at g = if we take the 
incommensurate term in Eq. to be of the form 5cos{Qn), where Q — Qo + q and Qq/it is equal to 
some simple rational number different from 1 . Let us first see this for the case of noninteracting fermions. 
Imagine filling up the Fermi sea up to a Fermi energy Ep = —coskp where kp = the Fermi 

velocity is then v — sinfc^^. We define the continuum Dirac field as 

il){x) = e:xjp{ikpx) + ipLix) exp(— jfci?a:) , (48) 

After dropping terms like exjp{±ikpx) which vary rapidly, the Dirac equations take the form 
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i ( d/dt + vd/dx ) ipR + - exp{iqx — ikp) ipi^ = 
i ( d/dt ~ vd/dx ) tpL + exp{—iqx + ikp) tpR = 



(49) 



The solutions of this equation have the plane wave form 



/ TpR{x,t) \ 



exp 



( 'P/"! { '^R expiikx + iqx/2 — ikp/2) \ 
\ flL exp{ikx — iqx/2 + ikp /2) J 



(50) 



On substituting this in (Ha), we find that there is a single energy gap lying between and given by 



E± = ^{vq ± S 



(51) 



Thus the size of the gap, 6, is independent of q. (All these calculations assume that q and S are both 
small; the gap size may depend on g if g is not small). If S is held fixed and q is increased from zero, a 
state first appears at zero energy when q reaches the value S/v. It is instructive to express all this in the 
language of bosonization. Since K = 1, we can combine Eqs. (^0|) and (^) to show that the sine-Gordon 
Lagrangian takes the form 



C 



1 

2'KV 



{d(j)/dtf - {d<i>/dx) 



16v 



[1 — cos(2(/) + ga; — /cp) 



(52) 



Note that we have fixed the coefficient of cos(2(/)) in such a way that, for g = 0, the soliton mass, 
including the quantum correctional3, is exactly equal to the gap from zero energy, i.e., 5/2. On shifting 
by {qx — kp)/2, we get 



C 



1 

2TTV 



{d(j)/dtf ~ (d(j)/dx - q/2f 



16w 



[ 1 - Q.O^(24>) 



(53) 



The Hamiltonian therefore contains a boundary term 



= - 1^ [ <t>{^) - 0(-oo) ] (54) 



which is purely topological; it is equal to —qv/2 in a one-soliton state. We thus see that if q is increased 
from zero, the energy of the one-soliton state becomes equal to that of the ground state (whose soliton 
number is zero) when q reaches the value S/v. From (^ij), we observe that a soliton corresponds to a hole; 
thus the energy of a hole (i.e., the energy required to remove one fermion from the system) becomes zero 
at q = S/v. Conversely, the energy required to add a fermion (antisoliton) to Ihfi^system becomes zero 
at q = —5 /v. [Similar results appear in the context of incommensurate crystalajO'EI] . 

To illustrate these results, we show the two major bands and the big gap separating them near the wave 
number Qq — 27r/3 and the Fermi energy Ep = — cos(7r/3) = —1/2 in Fig. 5; we have taken 5 = 0.4. 
[We have chosen a bigger value of S here compared to Figs. 1 and 2 in order to discuss some additional 
gaps in Sec. V and Fig. 6 below; those small gaps would not have been visible in Fig. 5 if we had chosen 
values of S much less than 0.4]. We used the condition tr M{q,E) < 2 to find the energy levels, and we 
used the energy resolution dE = 10~^. Note that the upper and lower edges of the gap follow Eq. (|5l]), 
and the gap closes with respect to Ep at q/S = ±l/w, where 1/v = 2/^/3 — 1.1547.... It is remarkable 
how simple Fig. 5 is compared to Figs. 1 and 2. 

Thus, if Qo ^ TT, there is a critical and nonzero value of q at which the gap closes at zero energy. This 
is very different from the earlier case with Qo = ""j where there is always a band of energies lying around 
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zero energy, so that there is no gap at zero energy for any q ^ Q. (Our resuhs disagree with Ref. 8 which 
argues that there is a nonzero critical value of q for Qq — • 

Finally, we can show that there is a nonzero critical value of q where the gap closes if Q 7^ tt and the 
fermions are interacting. By following arguments similar to the ones above, we find that the continuum 
theory is described by a combination of the Lagrangian in (^) (for q = 0) and the boundary term in 
(^) (for q nonzero but small). The soliton gap for g = is given by /S.E ^ ,51/(2-^); hence the boundary 
term implies that the soliton (antisoliton) gap will become zero at g = ±qc where 

qc = • (55) 

V 

We thus get an expression for the critical value of q in terms of the q = Q gap and the two Luttinger 
parameters K and v. 



V. HIGHER ORDER CONTINUUM THEORIES FOR NONINTERACTING FERMIONS 

WITH INCOMMENSURATE HOPPING 



In this section, we will consider the complete energy spectrum of noninteracting fermions with a weak 
incommensurate hopping. (This will complement the pictures of the complete spectrum of noninteracting 
fermions. with an incommensurate on-site potential given in the works of Hofstadter and of Sun and 
RalstonQ). In contrast to Sees. II and IV, we will use continuum theories which go up to higher than 
linear order in 5. This will enable us to look at some of the smaller gaps in the energy spectrum as can 
be seen in Fig. 7. 

The basic point is the following. Suppose that the incommensurate term Vg has Fourier components 
at wave numbers ±Q (where Q — Qq + q), and that we are interested in the vicinity of the two Fermi 
points at momenta ±kp. (We assume that < Qo < 27r and < kp < tt). Then a gap develops at the 
Fermi points if there exists an integer m (positive or negative) such that 

— kp + mQo — kp mod 27r . (56) 

To lowest order in the strength of the incommensuration S, the gap at the Fermi points is of order j'™', 
since Vs must act \m\ times to take us from the momentum state —kp to kp or vice versa. Now suppose 
that there is another integer n such that 

— kp + nQo = kp mod 27r . (57) 

We now have the possibility of interference between the actions of V^'™' and vj"', i.e. between terms in 
the continuum theory of the forms jl^le*""?^ and 5l"le"'^^. This interference gives rise to a large number 
of gaps and bands as we can show through a harmonic potential approximation. It is clear that conditions 
(^) and ( p7[ ) can hold simultaneously only if Qq/tt and kpjix are rational numbers. 

There are now two possibilities which we will call A and B. In case A, m = — n, which necessarily 
means that kp — 11/2 and Ep = 0. In that case, all the gaps are of the order j'™'. The discussion in Sec. 
II B and Figs. 1 and 2 is an example of case A, with Qq = tt, kp = tt/2 and m = —n = 1. In case B, 
\m\ < \n\. Here, there is one big gap of order jl™' and all the other gaps are of order j'"'. 

As an example of case B, let us consider the two small gaps originating from the top and bottom of 
the big gap shown in Fig. 5 at q/S = 0. We will only discuss the gap at the top here. An enlargement 
of that region is shown in Fig. 6 where it is clear that there are actually several small bands and gaps 
emerging from the top of the big gap; only the largest of these small gaps is visible in Fig. 5. This region 
can be understood as follows. Following the notation in Eqs. ( p6[]57| ), we have Qo = 27r/3, kp = 7r/3, 
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m ~ 1 and n = —2. As discussed in Sec. IV, the effect of ( |56| ) is linear in S. We now have to consider 
the effect of (|5^) which is of second order in S. The fact that two applications of Vg can take us from the 
modes near —kp to the modes near kp has three different effects. 

(i) At g = 0, it shifts the position of the lower edge of the upper band from the value Eq=Q = Ep + = 
— 1/2 + (5/2 given in Eq. ( pi] ) with Ep = —1/2. This shift can actually be calculated exactly, not just to 
order 6^, by considering the three states with momenta — 7r/3, 7r/3 and tt and writing down the matrix 
elements of Hq + Vs between these states as a 3 x 3 matrix. One of the eigenvalues of this matrix gives 
the position of the lower edge of the upper band to be 

1 A 1 / N 1/2 

Eg=o ^ - + - - - [9 ^ 66 + 3d' ) . (58) 

(ii) At second-order in Vs, we can go from the a state near — 7r/3 to a state near tt and then back to the 
original state near — 7r/3; similarly, we can go from 7r/3 to tt and back to 7r/3. This changes the velocity 
V = dE/dk from v — V3/2 to 

(iii) Finally, we can define the continuum Dirac fields as in (p8[). After performing the transformation in 
(pG) and shifting x by an appropriate amount, we get the time- independent equations 

S 5 

[ E + ivd/dx ]aR + -[!-— exp(-i3gx) ] = , 

S 5 

[ E — ivd/dx ] flL + 7: [ 1 — 77: exp{i3qx) ] an = , (60) 
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where 



E ^ E - ^ - E,^o + P , 

/3=^(1-^). (61) 

Thus we get a periodic potential as in Sec. II B. Following the procedure in that section, we first study 
the harmonic oscillator energies near the bottom of one of the wells, say, at a; = 0; then we will briefly 
comment on the broadening of those energies due to tunneling between different wells. On expanding 
exp{zLi3qx) up to order g^, wc find that the fields a± — aj^zLuL satisfy the harmonic oscillator equations 

[ - I d'/dx' + i u;'x' + ] a± = , (62) 



where the harmonic frequency is 



^ ^ (63) 

4 V 



The harmonic approximation is good if the root mean square deviation of the particle {x')-^/' ^ 1/v^ is 
much smaller than 27r/3g which is the periodicity of the potential in (|60|). We therefore assume that 

V \q\ « <53/2 . (64) 
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We then find the energy levels to be 



vq 



1 , vqS 



En = Eq^„ + ^ ^ f3 + (3' + 2v'uo{n + -) + , (65) 



2' 



1/2 



where n > 0. The term vqd'^/8 inside the square root comes from perturbatively including the first 
anharmonic correction. The eight solid lines in Fig. 6 show the energies in ( |65|) for n = 0,1, 2, 3, and for 
positive and negative values of q. The agreement is good for small values of q. We can now use the WKB 
approximation to study the band widths due to tunneling between wells; we find that the widths go as 

4 (5^/2 

Ai?„ ~ exp( - g ^ ) • (66) 



If 6 and q/5 are both small, then ( pq) is much larger than the factor of exjp{—2S/q) obtained in Eq. (|21 
Qualitatively, this explains why the bands in Fig. 6 broaden much more rapidly than the ones in Fig. 2 
as q increases from zero. 

Fig. 7 shows the complete picture of bands and gaps for all values of Q with S = 0.4. We see that the 
gaps come in two classes. As discussed above, a gap arising from case B consists of a small gap (which, 
at higher resolution, is actually a group of bands and gaps) emerging from one point on one edge of a 
big gap. On the other hand, a gap arising from case A looks like a cross in which four groups of gaps 
meet at a point; this can be seen in Fig. 7 at £' = and Q/n = 1/m where m = 1,2, 3, .... The simplest 
example of this corresponds to m = 1 which was discussed in Sec. II B. 

Let us now consider higher values of m. In general, we can go from the states near — 7r/2 to states near 
7r/2 in m steps, either through the intermediate states with increasing momenta at kj = — 7r/2 + jir/m 
with j = 1, 2, TO, or through the intermediate states with decreasing momenta at kj — — 7r/2 — jir/m 
with j = 1,2,..., TO. In each case, we get a product of to matrix elements from (||) divided by to — 1 
energy factors like Ej — Ep where Ep = 0. To be explicit, this product is 

nro'sin[(p+i)V"] <5™ 
= (Snr '-,n-i . = (67) 



n;ri sin[p^/TO] m2- 



(The second equality follows from the identity 

n— 1 n—l n 1 

n[2sin(p7r/n)] = lim TT (z - e'^^^/") = lim ^— = n ). (68) 

p—1 p—1 

After adding the contributions from the two sets of intermediate states and shifting x appropriately, we 
get a term in the continuum Hamiltonian of the form 

2F,n ( tplfi^L + i^li^R ) cos{mqx) . (69) 



This gives rise to a Dirac equation whose form is identical to that in Eq. ( |12[ ) except that q and S in that 
equation are now replaced by mq and 

(5™ 

2F™ = — (70) 

TO 2™-i ^ ' 

respectively. We therefore get the same pattern of bands and gaps as in Figs. 1 and 2; we recall here that 
those figures only show one-fourth of the region of interest since they are restricted to positive values of 
q and E. The complete regions which resemble crosses are visible in Fig. 7 along the line E = at the 
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values Q/tt = l/m. The change of scale from S to 2F,n in Eq. ( [70|) explains why the sizes of the crosses 
decrease rapidly as m increases. 

This concludes our discussion of the main features of the complete energy spectrum with incommen- 
surate hopping. Due to the close relationship of our work with the anisotropic AH problem, it is not 
surprising that our Fig. 7 resembles the complete gnergy spectrum with an incommensurate on-site 
potential as given in, say, Fig. 7 of Sun and RalstonQ. (Their variables $/$o and A correspond to our 
Q/{2n) and 1/d respectively). For instance, their figure also shows a series of crosses along the line E — 0. 
Due to the differences in the matrix elements in (^) versus the ones in (|37|), the scales of the crosses go as 

(5™ 

2F^ = — (71) 

m 

for the case of on-site incommensuration. This decreases more slowly than ( [70| ) with increasing m which 
explains why the crosses in Sun and Ralston are bigger than the ones in our Fig. 7; they also have larger 
values of S equal to 1 and 2/3 compared to our S = 0.4. Similarly, their figure shows the smaller gaps on 
the edges of the larger gaps as explained by our case B. 



VI. OUTLOOK 

The continuum method discussed in this paper is robust in the sense that it can be applied even if 
we add weak next-nearest-neighbor hoppings or the incommensuration is given by the sum of several 
harmonics. This is an advantage over_some other analytical methods. For instance, the Bethe ansatz 
discussed by Wiegmann and ZabrodinQ only works in somewhat special cases; in those cases, however, 
the Bethe ansatz works for all values of 6. A limitation of the continuum method is that it is accurate 
only if the incommensuration is weak. 

An interesting problem for future study may therefore be to examine what happens if the incommen- 
suration in the hopping is strong, i.e., if S is comparable to or larger than 1. [Both the perturbative and 
the continuum Dirac theory (or bosonization) approaches would then fail. Further, for strong incommen- 
suration, it does make a difference whether the incommensurate term is in the hopping or on-site; the 
latter case has been studied much more extensively in the literature] . The pattern of bands and gaps is 
expectedjiaJpg much more complicated in that case, perhaps with a devil's staircase or a point spectrum 
structureEfDU'Ll. The nature of the wave functions, namely, whether they remain extended or become 
localized, would also be of interest. 
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Figure Captions 



1. The bands as a function of the energy E and the wave number q (near dimerization), both in units of 
S, taking S = 0.05. The numbers 1, 3 and 5 labehng the three biggest gaps are explained in the text. 

2. A finer view of the bands as a function of E/6 and q/5, for S = 0.05. A comparison with the 
second-order perturbation results is indicated by the six solid lines. 

3. Semi-log plot of the widths AEn/S of the lowest four bands n = 0, 1,2,3. The points indicate the 
values obtained numerically for d = 0.05, while the solid lines (labeled by n) show the WKB expressions 
in Eq. (|l|). 

4. The density of states p as a function of E/S near dimerization. The solid and dashed lines show the 
behavior for q —^ and q — respectively. 

5. The bands as a function of {E — Ef)/5 and q/6, for Qo ~ 27r/3 and S — 0.4, with resolution dE ~ 10^^. 

6. Details of the bands as a function of {E — Ef)/S and q/S, near the top of the Qo ~ 27r/3 gap for 
S = 0.4, with resolution dE = 10"**. 

7. The bands as a function of E and Q/{2n) for 5 — 0.4, with resolution dE = 10""'. 
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FIG. 1. 




FIG. 2. 
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